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Abstract 

The generating function that records the sizes of directed circuit par- 
titions of a connected 2-in, 2-out digraph D can be determined from the 
interlacement graph of D with respect to a directed Euler circuit; the 
same is true of the generating functions for other kinds of circuit parti- 
tions. The interlace polynomials of Arratia, Bollobas and Sorkin [J. Com- 
bin. Theory Ser. B 92 (2004) 199-233; Combinatorica 24 (2004) 567-584] 
extend the corresponding functions from interlacement graphs to arbi- 
trary graphs. We introduce a multivariate interlace polynomial that is an 
analogous extension of a multivariate generating function for undirected 
circuit partitions of undirected 4-regular graphs. The multivariate poly- 
nomial incorporates several different interlace polynomials that have been 
studied by different authors, and its properties include invariance under a 
refined version of local complementation and a simple recursive definition. 

Keywords, circuit partition, interlace polynomial, isotropic system, 
local complementation, pivot, split graph 
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1 Introduction 

A graph G consists of a finite set V(G) of vertices, and a finite set E(G) of edges: 
each element of E(G) is incident on one or two vertices. An edge incident on only 
one vertex is a loop. It is often convenient to think of an edge as consisting of two 
distinct half-edges, each of which is incident on precisely one vertex. An edge 
is directed by specifying that one half-edge is initial and the other is terminal; 
as the half-edges are distinct, every edge can be directed in two different ways. 
The degree of a vertex is the number of incident half-edges, and a k -regular 
graph is one whose vertices all have degree k. Edges incident on precisely the 
same vertices are parallel, and a graph with no loops and no parallels is simple. 
A circuit in a graph is a sequence v\, hi, h' 2 , V2, Vk, hk, h' k+l , Ufc+i = v i 
such that for each i, hi and h[ are half-edges incident on Vi, and hi and h' i+1 arc 
half-edges of a single edge e^. A vertex may appear repeatedly on a circuit, but 
an edge may not appear more than once. If it happens that for every i, ej is a 
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Figure 1: C * v is obtained by reversing either of the two v-to-v walks within 
the incident circuit of C . 



directed edge with initial half-edge hi, then the circuit is directed] in general a 
directed graph may contain both directed circuits and undirected circuits. An 
Eulerian graph is a graph that possesses at least one Eulerian circuit, i.e., a 
circuit which includes every edge. 

This paper concerns a family of graph invariants, the interlace polynomials. 
We use the term in a generic sense, to include also some polynomials that were 
introduced under other names. All of these polynomials are motivated by the 
circuit theory of 4-regular graphs. 

Four cornerstones of this theory were laid in the 1960s and 1970s. A con- 
nected 4-regular graph is Eulerian, of course. More generally, an arbitrary 
4-regular graph has Euler systems, each of which contains one Euler circuit 
for each connected component of the graph. Kotzig 47 introduced the Hi- 
transformations: if C is an Euler system of a 4-regular graph F and v £ V(F) 
then the K-transform C * v is the Euler system obtained from C by reversing one 
of the two v-to-v walks within the circuit of C incident on v. Kotzig's theorem 
is the first of the four cornerstones; it tells us that all the Euler systems of F 
can be obtained from any one using ^-transformations. 

Although our discussion is focused on 4-regular graphs, we should certainly 
mention that Kotzig's theorem extends to arbitrary Eulerian graphs; see Fleis- 
chner's books [311 152"] for an account of the general theory. 

The second cornerstone of the circuit theory of 4-regular graphs is the in- 
terlacement graph I(F, C) of a 4-regular graph with respect to an Euler system 
C. I(F, C) has the same vertices as F, and v ^ w <E V(F) are adjacent in 
I(F, C) if and only if they are interlaced with respect to C , i.e., they appear in 
the order v...w...v...w on one of the circuits of C. The graphs that arise as in- 
terlacement graphs are called circle graphs. (This definition is usually restricted 
to Euler circuits of connected 4-regular graphs, but the restriction would be in- 
convenient here because there are natural ways to recursively simplify 4-regular 
graphs, which sometimes disconnect them.) This construction was discussed by 
Bouchet pj] and Read and Rosenstiehl [59] , who observed that the relationship 
between T(F, C) and I(F, C *v) is described by simple local complementation at 
v: iiv=/=x^=y^=v and x, y are both neighbors of v in F then they are adjacent 
in X(F,C * v) if and only if they are not adjacent in I(F,C). Later, Bouchet 
introduced isotropic systems to study circle graphs and the equivalence relation 
on arbitrary graphs generated by simple local complementations |T3l [14] [15] [17] . 

By the way, we use the term simple local complementation to distinguish 
this operation from the one that Arratia, Bollobas and Sorkin called local com- 
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Figure 2: D\ (top left) and D2 (bottom left) are nonisomorphic directed graphs 
whose undirected versions are isomorphic. They yield interlacement graphs that 
are equivalent under local complementations but not under pivots. 



plementation in [2] [3] [4] ; that operation also includes loop-toggling at neighbors 
of v. 

If C is an Euler system of F, then F is made into a 2-in, 2-out digraph D 
by choosing either of the two orientations for each circuit of C, and directing 
the edges of F accordingly. If v ^ w are neighbors in T(F, C) then the iterated 
K-transform C * v * vo * v is also a directed Euler system for D, obtained by 
interchanging the two v-to-w walks within the incident circuit of C . Following 
[21 [3], we refer to the operation C^C*v*w*vas transposition; the induced 
operation on interlacement graphs is the pivot G >-> G vw . Kotzig [47], Pevzner 
[55] and Ukkonen [72] proved that transpositions suffice to obtain all the directed 
Euler systems for a 2-in, 2-out digraph from any one; a more general form of 
this theorem was proven by Fleischner, Sabidussi, and Wenger [33 . 

As examples of these notions, consider the 2-in, 2-out digraphs D\ and D2 of 
Figure [2] They are small enough that each has only one directed Euler circuit, 
up to automorphism. Two Euler circuits are indicated in the figure; to trace an 
Euler circuit follow the directed edges, making sure to maintain the same dash 
pattern when traversing a vertex. (The dash pattern may be changed while 
traversing an edge.) The corresponding interlacement graphs are indicated in 
the figure's third column. Pivoting on an edge in the lower interlacement graph 
produces an isomorphic replica, with a different degree-2 vertex; pivoting on an 
edge in the upper interlacement graph has no effect at all. The fact that the two 
interlacement graphs are not equivalent under pivots reflects the fact that D\ 
and D2 are not isomorphic. On the other hand, simple local complementation at 
the single degree-2 vertex of the lower interlacement graph produces the upper 
interlacement graph, reflecting the fact that the undirected versions of D\ and 
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Figure 3: An orientation of the circuit of C incident at v is indicated at left. 
Each circuit partition of F involves one of the three pictured transitions at v. 
Both orientations of the incident circuit of C result in the same labeling of these 
three transitions. 



D 2 are isomorphic. 

Let F be a 4-regular graph with c(F) connected components, and let C be an 
Euler system of F. A circuit partition or Eulerian partition of F is a partition 
of E(F) into edge-disjoint circuits. Such a partition is determined by choosing, 
at each vertex of F, one of the three transitions (pairings of the incident half- 
edges): the transition that appears in the incident circuit of C, which we label 
<fi, for "follow" ; the other transition consistent with the edge-directions given by 
the incident circuit of C, which we label x, for "cross"; or the transition that 
is inconsistent with these edge-directions, which we label ip. See Figure [3] (We 
should mention that we use the terminology of Ellis-Monaghan and Sarmiento 
[29] and Jaeger [41] . in which a transition at v specifies both pairings of incident 
half-edges that might appear in a circuit partition. Other authors, including 
Bouchet and Kotzig, use "transition" in a slightly different way, to refer to a 
single pairing of half-edges, and require a separate matching-up of the pairings.) 
If n = \V(F)\ then F has 3" circuit partitions, given by choosing one of the 
three transitions at each vertex. A 2-in, 2-out digraph has 2™ directed circuit 
partitions. 

The third cornerstone of the circuit theory of 4-regular graphs is the idea 
of defining a polynomial invariant of a 2-in, 2-out digraph (or 4-regular graph) 
by using some form of the generating function ^ a;' p ' that records the sizes of 
(un)directed circuit partitions. This idea was introduced by Las Vergnas [51] . 
who observed that a polynomial defined recursively by Martin [54] is essentially 
equivalent to the generating function. Las Vergnas also extended the idea to 
general Eulerian graphs, and to 4-regular graphs that arise as medial graphs 
imbedded in the projective plane and the torus (49] [50]. In particular, if a 4- 
regular graph F is imbedded in the plane then its complementary regions can be 
colored checkerboard fashion, yielding a pair of dual graphs with F as medial, 
and the Tutte polynomial of either of the two dual graphs yields the directed 
circuit partition generating function for a certain directed version of F; this 
theorem foreshadowed the famous connection between the Tutte polynomial 
and the Jones polynomial of knot theory [42] [44] [63] . We will not focus any 
further attention on imbedded graphs in this paper; we refer the interested 
reader to Ellis-Monaghan and Moffatt [28] for a thorough discussion including 
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Figure 4: The two digraphs obtained by directed detachment at v. 



recent results. 

Martin observed that the circuit partition generating functions can be de- 
scribed recursively. Suppose F is a 4-regular graph with an Euler system G, 
D is the 2-in, 2-out digraph corresponding to a choice of orientations for the 
circuits of G, and v is unlooped in F. The directed circuit partitions of D fall 
into two classes, those that follow G through v and those that involve the x 
transition at v. These two classes correspond to directed circuit partitions of 
the two digraphs and D x obtained by directed detachment at v, illustrated 
in Figure HJ (The term detachment was coined by Nash- Williams; see [57] for 
instance.) 

Similarly, the circuit partitions of F fall into three classes according to the 
transitions at v, and the three classes correspond to circuit partitions of three 
graphs F$, F x and F^ obtained by detachment at v. 

The fourth cornerstone involves an equality due to Cohn and Lempel [22] . In 
its original form, the equality relates the number of cycles in a permutation to 
the nullity of an associated skew-symmetric matrix over GF(2). An equivalent 
form of the equality relates the number of circuits in a directed circuit partition 
of a connected 2-in, 2-out digraph D to the nullity of the adjacency matrix of 
an associated subgraph of an interlacement graph I(D,C). It is remarkable 
that versions of this useful equality have been discovered and rediscovered by 
combinatorialists and topologists so many times [6l f8] fl2 l rT9l l40l [43l l46l f52l f53l 
155] l56l 160] |6TJ |62] [73] . The Cohn-Lempel equality extends to a circuit-nullity 
formula for undirected circuit partitions in undirected 4-regular graphs. We 
simply state the formula here, and refer to [521 1Z0] for more detailed accounts. 
Let G — 1(F, G), let P be a circuit partition of F, and let Gp = Ip(F, C) be 
the graph obtained from G by removing each vertex at which P involves the <p 
transition used by G, and attaching a loop at each vertex where P involves the 
ip transition. Then the circuit-nullity formula states that 

\P\-c(F) = u(G P ), 

where v{Gp) denotes the GF(2)-nullity of the adjacency matrix of Gp, i.e., the 
difference between |V(Grp>)| and the GF(2)-rank. (The adjacency matrix of Gp 
is the V(Gp) x V(Gp) matrix over GF(2) in which a diagonal entry is nonzero 
if and only if the corresponding vertex is looped, and an off-diagonal entry is 
nonzero if and only if the two corresponding vertices are adjacent.) The original 
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Cohn-Lcmpcl equality is essentially the special case in which c(F) — 1 and no 
tp transition appears in P. 

Some examples appear in Figure [3] In the top row we see three circuit 
partitions Pi, P2, P3 of the undirected version of the graphs D\ and D% of Figure 
[5J the transition labels indicate the relationships between these circuit partitions 
and the Euler circuit of D\ indicated in Figure^ The corresponding graphs Gp i 
appear in the second row. In the third row of Figure [5] we see the same three 
circuit partitions, now with transition labels that indicate their relationships 
with the Euler circuit of D 2 indicated in Figure [5J The corresponding graphs 
Gp i appear in the fourth row. The circuit- nullity formula is satisfied because 
both sets of graphs Gp i satisfy u(Gp 1 ) =2 = 3 — 1, v(Gp 2 ) = 1 = 2 — 1 and 
is(Gp 3 ) =0 = 1-1. 

If F is a 4-regular graph with an Euler system C, then the circuit-nullity 
formula provides a bijective equivalence between these two information-sets. 

1. the circuit partitions of F; each circuit partition is determined by choosing 
the 4>, x or ij) transition at each vertex, and is listed with the number of 
circuits it contains 

2. the looped full subgraphs of I(F, C); each looped full subgraph is obtained 
by choosing to remove, retain or loop each vertex, and is listed with the 
GP(2)-nullity of its adjacency matrix 

This bijective equivalence tells us that the generating functions that record 
the sizes of (un) directed circuit partitions in F - or equivalently, the Martin 
polynomials - can be reformulated as generating functions that record the bi- 
nary nullities of adjacency matrices of (looped) full subgraphs of C). Ob- 
serve that there is no reason the nullity-based reformulations of these generating 
functions should be restricted to interlacement graphs; the definitions extend 
unchanged to arbitrary simple graphs. The graph polynomial that extends the 
directed Martin polynomial, the (vertex-nullity) interlace polynomial Qn{G), 
was introduced by Arratia, Bollobas and Sorkin [2J [3] . Their original definition 
extended the recursive description of the Martin polynomial, rather than the 
nullity-based reformulation; they derived the nullity-based form in g], where 
they also introduced a two- variable interlace polynomial q(G). Subsequently, 
Aigner and van der Hoist [TJ defined another interlace polynomial Q(G), which 
extends the reformulated version of the undirected Martin polynomial. More 
recently, Courcelle [26] introduced a multivariate interlace polynomial C(G), 
which extends a logically equivalent form of the information described in item 2 
above. These graph polynomials are related in various ways to the (un)restricted 
"Tutte-Martin polynomials" of isotropic systems studied by Bouchet [THl HB] . 

Our purpose in the present paper is to incorporate the </>, x, ip labels into 
the machinery outlined above. This is accomplished in two stages. In Section 
2, (j), x and ip are introduced into the circuit theory of 4-regular graphs as tran- 
sition labels. These transition labels are neither absolute nor arbitrary; they 
are defined with respect to a particular Euler system, and they are modified in 
particular ways when ^-transformations are applied to that Euler system. The 
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Figure 5: Three circuit partitions Pi, P2, P3 in the 4-regular graph of Figure [2 
with the transition labels and Gp graphs corresponding to the two Euler circuits 
shown there. 



7 



labeled versions of the Martin polynomials are generating functions that record 
the sizes of (un)directed circuit partitions, along with the corresponding tran- 
sition labels. (The V label is not needed for directed circuit partitions.) These 
labeled versions of the Martin polynomials are specializations of the transition 
polynomials discussed by Jaeger [UJ and Ellis- Monaghan and Sarmiento [25] ; 
the transition polynomials also incorporate transition labels, but the labels are 
arbitrary and consequently carry less information than </>, \ and tp. 

The equivalence between items 1 and 2 above tells us that the labeled ver- 
sions of the (un)directed Martin polynomials can be reformulated as labeled 
generating functions that record the GF(2)-nullities of adjacency matrices of 
(looped) full subgraphs of interlacement graphs. In Sections 3-5, these generat- 
ing functions are extended from labeled interlacement graphs to general labeled 
graphs, with vertex labels 0, x, V' now representing three ways to treat a vertex 
(remove it, retain it, or attach a loop) rather than three ways to choose a transi- 
tion. The result is to unify all the polynomial invariants of graphs and isotropic 
systems discussed in the paragraph before last in a single multivariate graph 
polynomial, the labeled interlace polynomial Q\(G), whose properties include a 
three-term recursive definition and invariance under labeled local complemen- 
tations. Setting -0 = in Q\{G) we obtain the 2-label interlace polynomial 
q\(G), which extends the nullity-based reformulation of the directed Martin 
polynomial. It satisfies a two-term recursion and is invariant under labeled 
pivots. 

In Section 6 we detail the relationships between Q\ and the several kinds 
of interlace polynomials that have been studied since the original definition of 
Arratia, Bollobas and Sorkin [2] [3] 0] . In Section 7 we present formulas for the 
labeled interlace polynomials of graphs with split decompositions, and discuss 
their computational significance. 

Acknowledgments. Before proceeding to present these notions in detail, 
we should thank D. P. Ilyutko, V. O. Manturov and L. Zulli. The idea of us- 
ing label-switching local complementations was inspired by many conversations 
with them while we studied the use of interlacement to describe the Jones poly- 
nomial and Kauffman bracket of a link diagram [37 ] 155 ] 155 ] 155 ] 157 ] IBS] [71] . 
Preliminary drafts of the paper were improved by the kind advice of B. Bol- 
lobas, B. Courcelle, J. A. Ellis-Monaghan, C. Hoffmann, M. Las Vergnas and 
an anonymous reader. 

2 Transition labels and circuit partitions 

Definition 1 An Euler system C in a ^-regular graph F is labeled by giving a 
trio of functions 4>c, XC ipe mapping V(F) into some commutative ring R. 

When we want to specify the ring in question, we refer to C as R-labeled. We 
will see in Section 6 that it is useful to consider a variety of labeling strategies, 
in a variety of rings. However, there is an especially natural way to implement 
Definition[T] Suppose R is a polynomial ring with 3- |^(-F)| indeterminates, one 
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indeterminate corresponding to each transition at each vertex of F. For each 
v G V{F), define the images of v under <pc, Xc, and ipe i n accordance with 
Figure [31 Then the label function 4>c actually specifies the Euler system C. 

Recall that a K-transformation is applied by reversing one of the two v-to-v 
walks in an Euler system, as depicted in Figure [TJ Clearly this reversal affects 
some transition labels: at v itself, the (j> and tp labels are interchanged; and at 
a vertex w that appears precisely once on the reversed v-to-v walk, the x an d 
V> labels are interchanged. 

Definition 2 Let v be a vertex of a ^-regular graph F , and let C be a labeled 
Euler system of F. The labeled /c-transform C * v is obtained by reversing one 
of the two v-to-v walks within the circuit of C incident on v, and making the 
following label changes: <ftc*v(v) — ijjc{ v ), ipc*v(v) — <Pc(v), and for each w ^ v 
that neighbors v in X(F,C), xc*v(w) — ipc(w) an d ipc*v(w) = Xc(w). 

Definition 3 Two labeled Euler systems of a ^-regular graph are K-equivalent if 
and only if one can be obtained from the other through labeled ^-transformations. 

Kotzig's theorem [47] tells us that ^-transformations can be used to obtain 
all the Euler systems of F from any one; it follows that if F has Euler systems 
C\ and C2 then each labeled version of G\ is /{-equivalent to a unique labeled 
version of C2. 

Definition 4 Let C be an R-labeled Euler system of F , and suppose y £ R. For 
each circuit partition P of F, let (j){P,C), x(P, C) and ip(P,C) denote the sets 
of vertices of F where P involves the transition labeled <f>, x or ip (respectively) 
with respect to C . The i?-labeled circuit partition generating function of F with 
respect to C is 



k(f,c)= ]t n n xc(v)\\ n ^ 



y \p\-cm t 



where V{F) is the set of circuit partitions of F and c{F) is the number of 
connected components of F. 

Different systems of labels yield generating functions with different lev- 
els of detail. If the labels take the natural values in the polynomial ring 
^[{y} U {</>c(v),Xc(v),i>c(v) I v G V(F)}} with I + 3 • \V(F)\ independent 
indeterminates, then tt(F, C) is essentially a table that lists, for every circuit 
partition P, \P\ — c(F) along with the relationship between C and P at every 
vertex. Similarly, the ip = specialization of this polynomial is essentially a 
generating function for directed circuit partitions of a 2-in, 2-out digraph D ob- 
tained by directing the edges of F according to orientations of the circuits of C. 
The theory of this polynomial is outlined in Section 5. If <f>, x, tp = 1 G Z[y] then 
tt(F, C) is the ordinary (unlabeled) generating function for circuit partitions of 
F, and if ip = and </>, x = f G Z[y] then ir(F,C) is the ordinary generating 
function for circuit partitions of D. 
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Figure 6: All Euler circuits have the same x transitions at the central vertices. 



Proposition 5 Let C\ and C% be labeled Euler systems in F. If C\ and Ci are 

n-equivalent, then ir(F, C\) = tt(F, C2). 

Proof. Labeled ^-transformations preserve tt term by term, because the 
contribution of each P G V(F) is unchanged. ■ 

The labeled circuit partition generating function ir(F, C) incorporates more 
precise information about the structure of a 4-regular graph F than other 
transition-based polynomials that have appeared in the literature. For instance, 
the "Tutte-Martin polynomials" discussed by Bouchet [16] involve "coding" each 
vertex of F by choosing an arbitrary labeling of the transitions. (Indeed Propo- 
sition (5.2) of 13! states that the isotropic systems associated to two coded 
4-regular graphs are isomorphic if and only if they can be obtained from dif- 
ferently coded versions of the same 4-regular graph.) Jaeger's transition poly- 
nomial |41j and the generalized transition polynomial of Ellis-Monaghan and 
Sarmiento 30J also involve arbitrary transition labels. In contrast, the label 
functions <pc , Xc and ipc are not arbitrary: they are associated in special ways 
with the positioning of C within F, and they are handled in special ways by 
labeled ^-transformations. 

In particular, it is not generally possible to simply transpose two labels at one 
vertex using labeled ^-transformations. For example, suppose C is any labeled 
Euler circuit of the graph F in Figure El with 21 independent indeterminates 
serving as transition labels (three for each vertex). Then ir(F,C) determines 
the x labels at the two central vertices of the graph: they are the only transi- 
tion labels that appear only in terms divisible by y, because the corresponding 
transitions are the only ones that do not appear in any Euler circuit. Similarly, 
if F\ is the undirected version of the graph D\ of Figure El and C\ is the Euler 
circuit of D\ indicated in Figure El then the ip labels in F\ are distinguished by 
the fact that they are the only ones that appear in a term of n{F\, C\) divisible 
by y 3 - 

The labeled circuit partition generating function n(F, C) satisfies a labeled 
version of the detachment-based recursion mentioned in the introduction. Sup- 
pose C is a labeled Euler system of F and v is an unlooped vertex of F. (We 
leave the consideration of looped vertices to the reader.) Then as discussed in 
the introduction, there are three associated 4-regular graphs obtained by de- 
tachment at V, denoted F$, F x and F^ according to the transitions that define 
the detachments. As illustrated in Figure [3 C yields labeled Euler systems in 
all three detachments. F^ has a labeled Euler system whose circuits simply 
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Figure 7: The graph of Figure |6] appears on the left, with an Euler circuit 
indicated by dashes. The top row depicts the three detachments at a central 
vertex and the bottom row depicts the three detachments at another vertex. 



follow the circuits of C, omitting v. has a labeled Euler system Cy,, ob- 
tained in the same way from the labeled K-transform C * v. The situation in 
F x is more complicated, as there are two distinct cases. If v is not interlaced 
with any other vertex with respect to C, then c(F x ) — c(F) + 1 and F x has a 
labeled Euler system C x obtained by separating the two v-to-v circuits within 
the incident circuit of C . On the other hand if v is interlaced with a vertex w, 
then c(F x ) = c(F) and a labeled Euler system C x of F x is obtained by following 
the circuits of C * w * v * w, omitting v. 

Theorem 6 Suppose v is an unlooped vertex of F. Then 

n(F, C) = 4>c{v) ■ 7r(F , C ) + X c(v) ■ tt(F x , C x ) + ^ c {v) • 7r(F^, C^). 

Proof. The formula is justified in the natural way, by classifying the circuit 
partitions of F according to the transitions at v. The notation (F X ,C X ) is 
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ambiguous - it does not tell us whether the first case or the second case holds 
in F, and in the second case it does not reflect the fact that different choices of 
w will yield different labeled Eulcr systems in F x - but this ambiguity does no 
harm because the formula holds for every choice of C x . ■ 

3 Vertex labels and partitions 

As discussed in the introduction, the circuit-nullity formula allows us to reformu- 
late the theory described in Section 2: rather than thinking of circuit partitions 
in a 4-regular graph F, we may think of graphs obtained from an interlacement 
graph X(F, C) by removing some vertices and looping some vertices. This re- 
formulated version of the theory extends directly from interlacement graphs to 
arbitrary simple graphs. 

Definition 7 A simple graph G is (R-)\abeled by giving a trio of functions 
4>G,XGj' i pG mapping V(G) into some commutative ring R. 

Once again, the most natural ring to use as R is a polynomial ring with 
three independent indeterminates <Pg(v), Xg{v) and i/)g(v) for each v S V(G). 

Definition 8 Let G be a labeled simple graph with a vertex v. The labeled 
local complement G\ is the labeled graph obtained from G by toggling adjacencies 
between distinct neighbors ofv and making the following label changes: 4>g^{ v ) = 
4>g{v), tJjg v x { v ) — 0g(^); and for each w ^ v that neighbors v in G, xg«(u>) = 
ipG(w) and ip G i(w) = Xg («->)• 

Recall Definition [2j if C is a labeled Eulcr system in a 4-regular graph 
F, then for v € V(F) the labeled K-transform C * v is obtained by reversing 
one of the two v-to-v walks within the circuit of C incident on v, and making 
appropriate adjustments to the label functions. The effect of reversing a v-to-v 
walk is to toggle interlacements between vertices that appear precisely once on 
that walk, i.e., vertices that are interlaced with v. We conclude the following. 

Theorem 9 Let F be a ^-regular graph with a labeled Euler system C , and 
let T(F, C) be the corresponding labeled interlacement graph. If v G V(F) then 
l{F,C*v)=l(F,C)l. 

We now have three different kinds of local complementation: simple local 
complementation, for which we use no particular symbol; the local complemen- 
tation G v used by Arratia, Bollobas and Sorkin [H [3j H] , which combines simple 
local complementation with loop-toggling at neighbors of v; and labeled local 
complementation. The loop-toggling at neighbors of v in G v and the \ ^ "0 
exchange at neighbors of v in G x are essentially the same thing; what is special 
about labeled local complementation is the <j> ■<->■ ip exchange at v itself. As we 
show in Theorems [T3l and [TBI below, this detail allows us to formulate a simple 
invariance property for a rather complicated-seeming graph polynomial that de- 
termines all the different interlace polynomials studied in [TJ [2j [3j 01 [TBI US ES] . 
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It is the absence of this detail that creates the seeming lack of simple invariance 
properties in the original discussions of some of these graph polynomials. 

The reformulated version of a circuit partition is a certain kind of vertex 
partition. 

Definition 10 A labeled partition of a labeled simple graph G is a partition P 
of V(G) into three pairwise disjoint subsets, V(G) = <f>(P) U x(P) U ^(P)- The 
set of all labeled partitions of G is denoted V\{G). 

Note the different uses of the term "labeled" in Definitions [71 and ITUl a graph 
is labeled by specifying functions <j>Q, xg , i>G ■ V{G) — ¥ R, but a partition P 
of V(G) into three disjoint subsets is labeled by identifying one of the subsets 
as 4>(P), another as x(P), and the third as ip(P). 

Definition 11 If v G V(G) then the labeled local complement of P G V\(G) 
is the labeled partition P£ G P(G^) obtained from P by making the following 
changes: v G 4>(P\) if an d only ifvE ip(P), v G ip(P^) if and only if ' v G 4>(P)> 
and if w is a neighbor of v in G then w G x(P\) if an d only if w G ^(P), and 
w G tjj(P^) if and only if w € x(P)- 

The circuit-nullity formula suggests the following. 

Definition 12 If P is a labeled partition of a simple graph G then Gp denotes 
the graph obtained from G by removing every vertex in <j)(P) and attaching a 
loop at every vertex in tp(P). 

If P is a circuit partition of a 4-regular graph F with two Euler systems C\ 
and C2, then the circuit- nullity formula tells us that \P\ is related to the binary 
nullities of the adjacency matrices of the two Gp graphs obtained from X(F, C\) 
and T(F, C2). P itself does not change when we change Euler system, so these 
two adjacency matrices must have the same nullity. It may be a surprise that 
this invariance of nullity extends to arbitrary graphs, even though there is no 
fixed object that plays the role of P. 

Theorem 13 If v E V(G) then for every P e V\{G), v{G P ) = v((G x ) P v). 
Proof. If v G 4>(P)i the theorem states that 



and 



have the same GF(2)-nullity. Here the left-hand matrix is partitioned into sets 
of rows and columns corresponding respectively to the neighbors of v in x{P)i 
the neighbors of v in tp(P), and the vertices in x(P)U^(P) that are not neighbors 
of V] the first row and column of the right-hand matrix correspond to v. Bold 
numerals denote rows and columns with all entries the same, and an overbar 
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indicates the toggling of all entries in a matrix over GF{2). The nullity equality 
is verified by observing that adding the first row to every row in the first two 
sets of rows in the right-hand matrix yields 
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whose nullity is the same is that of the first matrix displayed above. 

If v € ip(P) the preceding argument is simply reversed. That is, we use row 
operations to show that 
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and 



have the same G.F(2)-nullity. 

If v € x{P) then the theorem states that 
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have the same GF(2)-nullity. This is verified by adding the first row to those 
in the second and third sets. ■ 

Here is the reformulated version of Definition |4j 



Definition 14 Let G be an R-labeled simple graph, and suppose y G R. 
labeled interlace polynomial of G is the sum 



The 



Qx(G) 




n ^(v) 



n * g ( w ) 

,V£x(P) 



Note that we call Q\ a polynomial even though y is an arbitrary element of 
R. This is a mere formality, as we could specify that y be an indeterminate and 
then obtain other instances of the definition through evaluation. 

Theorem 15 Let F be a ^-regular graph with a labeled Euler system C , and let 
I(F, C) be the corresponding labeled interlacement graph. Then 

n(F,C) = Q x (l(F,C)). 

As we mentioned in the introduction, the labeled interlace polynomial yields 
all the different kinds of interlace polynomials in the literature, by using different 
label values. This might suggest that the properties of Q\ would be more 
complicated than those of the other polynomials; instead the theory of Q\ turns 
out to be considerably simpler. 
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Theorem 16 Let G be a labeled simple graph with a vertex v. Then Q\(G) = 
Qx(GD- 

Proof. This follows immediately from Theorem Q2] and the fact that for 
every P G V X {G), 

n ^(«)) n xg(v)\ n ^ 
n ( n xgsm] ( n <m«) 

How can it be that Q\ is invariant under labeled local complementation, 
when choosing particular labels in Q\ yields the multivariable interlace poly- 
nomials of 4 and [26], which seem to have no invariance properties at all? 
The answer is given in Section 6 below: the choices of labels that yield these 
polynomials do not maintain the separation of 0, x an d ip- Losing this three- 
fold distinction makes it impossible to apply Definition [SI and consequently the 
properties of Q\ are obscured. 

Theorem 17 The labeled interlace polynomial of a labeled simple graph is re- 
cursively determined by these three properties. 

1. If G consists only of a single vertex v then 

Qx(G) = <Pg(v) + xg{v) ■ y + i>G{v)- 

2. If G\ and G2 are disjoint graphs then Q\(G\ U G2) = Q\(G\) ■ Q\{G2). 

3. Suppose v ^ w are adjacent vertices in a labeled simple graph G. Then 

Q X {G) = Mv)-Qx(G-v)+ XG (v)-Qx(((G^)l)-v)+ij G (v)-Q x ((Gl)-v). 

Proof. The first property is a special case of Definition [14J The second 
follows from the fact that a labeled partition P of G\ U G2 is simply the union 
of labeled partitions P%,P2 of G\ and G2 (respectively); consequently the adja- 
cency matrix of Gp is 

A x 
A 2 

where A\ and A2 are the adjacency matrices of (G\)p 1 and (G2)p 2 . 

The three summands in the recursive formula correspond to the natural 
partition V\{G) = S\ U S2 U ^3, with Si containing the partitions that have 
v G 4>(P)> S2 containing those that have v G x(-P), and S3 containing those that 
have v G i(>{P). 
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Consider the bijection / : Si ->■ V\{G - v) given by i>(f(P)) = ^(^), 
= X(-P) and <K/(P)) = 4>{P) ~ {v}. Definition QJ tells us that (G - 
v )f(P) ~ for every P e Si, so 4>g{ v ) • Q\(G — v) equals the sum of the 
contributions to Q\(G) of the partitions in S\. This justifies the first summand. 

According to Theorems H3] and [TBI if f G S3 then the contribution of P 
to Q\{G) equals the contribution of P£ to Q\(G^). Also, P G S3 if and only 
if v £ <I>{P\), so the argument just given for Si applies to {P% \ P G S3}. 
Consequently <j)G^(v) ■ Qa((G^) — v) equals the sum of the contributions to 
Q\{G V X ) of the partitions P£ with P G S3; as 4>gi(v) — iPg(v), it follows that 
V'g(v) • Q\((Gx) — v ) equals the sum of the contributions to Q\(G) of the 
partitions P G S3. This justifies the third summand. 

The second summand is justified in a similar way. If P G V\{G) then 
v G x(-P) if and only if « G if{P^ ), and this holds if and only if u G 0((P^)a); 
hence the argument given for Si applies to {(P\)x I P G S2}. ■ 

Unlike Theorem|51 Theorem[T7]provides a complete recursion. The difference 
is that looped vertices in 4-regular graphs, which are not covered in Theorem 
O give rise to isolated vertices in interlacement graphs, which are covered in 
Theorem [T7J 

The double local complementation in the second summand of part 3 of The- 
orem [17] may be replaced with a pivot; this version of the recursion is useful in 
proving Theorem 1261 below. 

Proposition 18 The formula of part 3 of Theorem ] 17\ may be rewritten as 
Q\(G) — 4>g( v ) • Q \{G — v) + xg{v) ■ Q\{{{{G\)x)\) ~ v ) +^g( v ) • Q\((G\) — v). 

Proof. (((G™)x) - U)™ = (((C?a)a)a) - V, so Theorem M tells us that the 
two formulas are the same. ■ 

4 Looped graphs 

Looped vertices play two very different roles in the theory discussed in the 
introduction. A 4-regular graph F may certainly have looped vertices, but in- 
terlacement graphs may not; they are simple by definition. Looped vertices 
reappear in the circuit-nullity formula, in association with ip transitions. Sim- 
ilarly, the definition of Q\(G) presumes that G is a labeled simple graph, but 
looped vertices play an important role because the graph Gp associated to a 
labeled partition P has a loop at each v G i(>(P). 

Observing that the difference between a looped vertex and an unlooped 
vertex in Gp is the difference between x(-P) and ip(P), we are led to a natural 
way to extend the theory of Section 3 to labeled, looped simple graphs. 

Definition 19 If G is a labeled simple graph with some looped vertices then its 
simplification is the labeled simple graph G s i mp obtained by interchanging xi v ) 
and ip(v) at each looped vertex, and then removing all loops. 
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The discussion of Section 3 is applied to a labeled, looped simple graph G 
indirectly, by using G S i mp as a stand-in for G. Note that with this approach 
looped graphs do not add anything new, so the difference between restricting 
the theory to simple graphs and extending the theory to looped graphs is es- 
sentially a matter of style, not substance. We choose to present the restricted 
theory because it is (appropriately) simpler. The extended theory requires more 
complicated statements of definitions and theorems - for instance, the extended 
version of Definition [5] would involve different label-swaps at looped vertices, 
and the extended version of Theorem [T7] would require an extra step to elimi- 
nate loops - and the complications seem unnecessary because all they amount 
to is the repeated application of Definition [T5] 

5 Labeled pivots and 2-in, 2-out digraphs 

Recall that if v is a vertex of a graph G then the open neighborhood of v in G 
is Ng(v) — {a G V(G) | a is adjacent to v in G}. 

Proposition 20 Let v and w be adjacent vertices in a labeled simple graph G. 
Then ((G%)l)f = ((GJQaOa = G ' 18 the lab e led simple graph obtained from G 
by making the following changes: 

1. 4> G '(v) = Xg(v), XG'(v) = <Pg{v) , (f>G'(w) = Xg(w) and xg>(w) = 4>g(w). 

2. Toggle the adjacency status of each pair of vertices a ^ b € V(G) — {«, w} 
such that a £ Ng(v), b € Nq(w) and at least one of a, b is not in Nq(v) n 
N G (w). 

3. Exchange the neighbors of v and w. 

Proof. This follows directly from Definition [8] ■ 

The operation G M> G' is the labeled pivot on the edge vw; we denote it 
G i ^ G\ w . Arratia, Bollobas and Sorkin [2j |3] noted that for unlabeled simple 
graphs, the result is the same up to isomorphism if the neighbor-exchange of step 
3 is replaced by a "label swap" in which the names of v and w are exchanged. 
An analogue of their observation holds here too: the result is the same up to 
isomorphism if steps 1 and 3 of Proposition [20] are replaced by the following 
swap of labels at v and w: 4*g'( v ) — Xg(w), XG'( v ) — 4>g{w) , ipG'(v) — iPg(w), 
<Pg>(w) = xg(v), XG'(w) = <Pg{v) and ipG>{w) = ^g(v). 

As mentioned in the introduction, the transposition operation ChC*w* 
v * w (where v and w are interlaced on C) plays the directed version of the 
role played by ^-transformation for undirected 4-regular graphs: if D is a 2- 
in, 2-out digraph then every directed Euler system of D can be obtained from 
any one through transpositions. Consequently pivots play the same role in the 
theory of directed interlacement as local complementations play in the theory 
of undirected interlacement. With this idea in mind, it is easy to formulate the 
following directed versions of the definitions and results discussed in Sections 2 
and 3. We leave the proofs to the reader. 
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Theorem 21 Let D be a 2-in, 2- out digraph with a directed Euler system C and 
let T(D, C) be the corresponding labeled interlacement graph. If v ^ w £ V(D) 
are neighbors in X(D, C) then X(D, C *w *v * w) — I(D, C) v x w . 

Definition 22 Let D be a 2-in, 2-out digraph, and let V(D) denote the set of 
directed circuit partitions of D. If C is a labeled, directed Euler system of D 
then for each P £ V{D), let 4>c(P) o,nd Xc{P) denote the sets of vertices of D 
where P involves the transition labeled <f> or \ (respectively) with respect to C. 
The labeled directed circuit partition generating function of D with respect to 
C is 

*{D,C)= J2 ( II ( II Xc(v) ) y\ p \-« D \ 

Pev(D) \ve<p(P) J \ve x (P) J 

where c(D) is the number of connected components of D. 

Definition 23 Let G be a labeled simple graph, and let V\o(G) = {P £ V\(G) 
| ip(P) = 0}. The 2-label interlace polynomial of G is 



Pev X o(G) \ve4>(P) J \ve x (P) J 

Theorem 24 Let G be a labeled simple graph with adjacent vertices v ^ w. 
Thenq x (G)=q x (Gl w ). 

Theorem 25 Let D be a 2-in, 2-out digraph with a directed Euler system C 
and let I(D, C) be the corresponding labeled interlacement graph. Then 

q x (l(D,C))=7r(D,C). 

Theorem 26 The 2-label interlace polynomial of a labeled simple graph is re- 
cursively determined by these three properties. 

1. If G consists only of a single vertex v then 

<7a(G) = 4> G {v) +xg(v) ■ y. 

2. If Gi and Gi are disjoint graphs then q\(G\ U G2) = q\(Gi) • q\(G2). 

3. Suppose v ^ w are adjacent vertices in a labeled simple graph G. Then 

q x (G) = 4> G (v) ■ qx{G -v) + XG (v) ■ qx((G v x w ) - v). 



The above results indicate that the entire theory of Section 2 can be re- 
stricted from 4-regular graphs to 2-in, 2-out digraphs, using directed Euler 
systems to describe directed circuit partitions via interlacement, and using la- 
beled transpositions and pivots rather than labeled ^-transformations and local 
complementations. This restriction is quite natural, but there is also a purely 
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algebraic way to restrict attention to directed circuit partitions, using inter- 
lacement with respect to arbitrary (undirected) Euler systems. The idea is to 
use the label functions to remove from consideration those transitions that are 
inconsistent with the edge-directions. 

Definition 27 Let D be a 2-in, 2-out digraph whose undirected version is F , 
and let C be a labeled Euler system of F. The label functions 4>c,XC,tpC are 
consistent with D if they have this property: For every vertex v £ V{F), the 
transition at v that is inconsistent with the edge- directions of D corresponds to 
a label value 4>c{v), xc(v), or ipc(v) that equals 0. 

The set of D-consistent labeled Euler systems of F is closed under labeled 
^-transformations, so D-consistent Euler systems provide a way to restrict at- 
tention to directed circuit partitions without any need to formally restrict the 
combinatorial machinery of Section 2. Of course the D-consistent versions of 
formulas like those of Definition [4] and Theorem [6] are simpler than they appear 
to be, because some of the label values are 0. 

5.1 T-compatible circuit partitions and Euler systems 

In this subsection we briefly discuss a notion introduced by Kotzig [4"T1 |4"8] and 
subsequently investigated by other researchers, including Fleischner, Sabidussi 
and Wenger [33], and Genest [341135]. 

Suppose F is 4-regular and T is a set that includes no more than one tran- 
sition at each vertex of F. The circuit partitions and Euler systems of F that 
avoid using the transitions from T are called T-compatible; Kotzig proved that 
T-compatible Euler systems exist for every T. For example, if D is a 2-in, 2-out 
digraph based on F then using T = {transitions inconsistent with the edge- 
directions of D} has the effect of restricting attention to directed circuits of D. 
Also, if we are given a (classical or virtual) link diagram, then the Kauffman 
states [SI [45] of the diagram are obtained by using T — {transitions corre- 
sponding to the strands of the link components incident at the crossings of the 
diagram}. 

Just as there are two ways to restrict the machinery of Section 2 to 2-in, 2-out 
digraphs, there are two ways to restrict the machinery to T-compatible circuits. 
The purely algebraic way involves the use of arbitrary Euler systems and T- 
compatible label functions, i.e., label functions with values corresponding to 
transitions in T. As the set of such labeled Euler systems is closed under labeled 
^-transformations, the full combinatorial machinery of Section 2 applies. 

The more thoroughgoing restriction involves considering only Euler systems 
that are themselves T-compatible, i.e., they do not involve any transition from 
T. The centerpiece of this restriction is the following T-compatible version of 
Kotzig's theorem, due to Fleischner, Sabidussi and Wenger [33] . 

Theorem 28 All the T-compatible Euler systems of F can be obtained from 
any one using these three types of operations: 
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1. n -transformations C t-> C * v, where T does not include a transition at v 

2. transpositions C^C*v*w*v, where T includes the tp transitions of C 
at v and w, and v and w are interlaced on C 

3. k -transformations C H- C*v, where T includes the \ transition of C at v 

Genest [34[ [35] codes the interlacement graph I(F, C) by coloring black the 
vertices that appear in part 2 of Theorem [551 an d coloring white the vertices 
of part 3. Then the ^-transformations and transpositions of Theorem 1281 yield 
local complementations and pivots that include color-swaps where appropriate. 
As discussed in the next section, the Arratia-Bollobas-Sorkin interlace polyno- 
mials of looped graphs are motivated by a different convention, involving the 
attachment of loops to the vertices of X(F, C) that appear in part 3 of Theorem 

on 

6 Interlace polynomials 

Arratia, Bollobas and Sorkin introduced the interlace polynomials after study- 
ing special properties of Euler circuits of 2-in, 2-out digraphs useful in analyzing 
DNA sequencing. The original interlace polynomial [2j [3] is a one- variable poly- 
nomial associated to a simple graph; following 0J, we denote its extension to 
looped graphs ^(G), and call it the vertex-nullity interlace polynomial. This 
polynomial was first defined recursively, using the local complementation and 
pivot operations. Using the recursive definition, Arratia, Bollobas and Sorkin 
proved that if G is the interlacement graph of a 2-in, 2-out digraph D then 
qN{G) is essentially the generating function that records the sizes of the parti- 
tions of E(D) into directed circuits. We refer to this fact as the fundamental 
interpretation of qN for circle graphs. 

In 0], Arratia, Bollobas and Sorkin showed that qn also has a non-recursive 
definition involving the nullities of matrices over the two-element field GF(2). 
For S C V(G) let G[S] denote the full subgraph of G induced by S, and let 
^(Gf.!?]) denote the nullity of the adjacency matrix of G[S] over GF(2). 

Definition 29 The vertex-nullity interlace polynomial of G is 

SCV(G) 

Considering Definition [121 we see that is obtained from Q\ by replacing y 
with y— 1, using (j),\ = \ and ip = at unlooped vertices, and using (f>, tp = 1 and 
X = at looped vertices. The basic theory of qN follows from this observation 
and the results of Section 3. (The <f>, x = 1 specialization of the polynomial 
q\ of Section 5 yields the restriction of q^ to simple graphs.) In particular, 
Theorem 1151 yields a fundamental interpretation of qN for looped circle graphs: 
Let F be a 4-regular graph with an Euler system G, and suppose G is obtained 
from I(F, C) by attaching loops at the vertices that appear in a certain subset 
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£ C V(G). For each circuit partition P e V{F), let (f>(P,C), x(P,C*) and 
V'(-P) C) denote the sets of vertices of P where P involves the transition labeled 
0, x 01 $ (respectively) with respect to C. Let V(F, £) denote the set of circuit 
partitions P with x(P C) l~l £ = and ip(P, C) C £. (That is, P(P, £) contains 
the T-compatible circuit partitions, where T includes the x transitions of C at 
vertices in £ and the ip transitions of C at vertices not in £.) Then 

9*r(G) = £ (y-lf^ F \ 

PGV(F,C) 

where c(P) is the number of connected components of F. 

Theorem [16] and Proposition [20] imply Remark 18 of 3 : if v and w are 
unlooped neighbors then qAr(G) = qN(G vw ). This equality does not extend to 
labeled pivots involving looped neighbors, because part 1 of Proposition [201 is 
not compatible with the special label values used to obtain from Q\. For 
the same reason, Theorem [12] does not yield a useful invariance property for qpj 
under local complementation at unlooped vertices. At looped vertices (j> = ip, 
though, so Theorem [TBI yields the equality gjv(G) = qN{G v ). Observe that 
according to the fundamental interpretation, if G is a looped circle graph then 
the equalities qN{G) = qN(G vw ) (for unlooped neighbors v, w) and Qn{G) = 
qN(G v ) (for looped v) follow immediately from Theorem |2"51 the T-compatible 
version of Kotzig's theorem due to Fleischner, Sabidussi and Wenger [33]. In 
general, these equalities indicate a connection between q^ and a well-known 
matrix operation, the principal pivot; detailed discussions are given by Brijdcr 
and Hoogeboom [50] [UJ and Glantz and Pelillo [55] . 

The two- variable version of the interlace polynomial was introduced in [?] : 

Definition 30 The interlace polynomial of a graph G is 

q(G)= {x-l)\ s \- v < G W\y-l) v <- G W\ 

SCV(G) 



x-l 



No fundamental interpretation was given for q in [1] , but rewriting Definition 
EOlas 

q(G)= ( X -^ S 

SCV(G) 

we see that q is obtained from the labeled interlace polynomial Q\ by using 
= 1, X = x — 1 and ip = at unlooped vertices, using 4> = 1, x = and 
ip = x — 1 at looped vertices, and replacing y with (y — l)/(x — 1). Consequently 
a fundamental interpretation of q for looped circle graphs follows immediately 
from Theorem 1151 Let F be a 4-regular graph with an Euler system C, and 
suppose G is obtained from T(F, C) by attaching loops at the vertices that 
appear in a certain subset £ C V(G). Then in the notation used above, 



E 

Pev(F,c) 



y~ 1 \ l rr _-i\\V{G)\-\4>{P,C)\ 
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Considering Definition [3] and Proposition[2ni we see why [J] does not mention 
any invariance properties of q under local complementations or pivots: because 
of label swaps involving 0, the exponent of a;— 1 associated to P € V{F, C) is not 
generally the same as the exponent associated to PV. Similarly, the two-term 
recursion for q that results from Theorem [T7] when we set one label to for each 
vertex does not appear in [3] because of the shifting-around of powers of x — 1 
under labeled local complementation. Some of these complications were handled 
in |64) by manipulating weights, but no motivation involving interlacement was 
provided there. 

Another interlace polynomial, denoted Q, was introduced by Aigner and 
van der Hoist in pQ. They showed that Q coincides with Bouchet's "Tutte- 
Martin polynomial" M [THJ [18]. In particular, Q has a simple fundamental 
interpretation for circle graphs: if G — T(F, C) then Q(G) is essentially the 
generating function that records the sizes of all the circuit partitions of F. The 
description in pQ makes it clear that Q is obtained from Q\ by using <f>, \i = !• 
As the labels are all the same, Theorem [TBI applies in this case; it tells us that 
Q is invariant under simple local complementation, as in Corollary 4 of [I] . 

Observe that the fundamental interpretation of q for circle graphs is quite 
different from those of qw and Q: only q has a fundamental interpretation as 
a generating function that records detailed information regarding the numbers 
of transitions of particular types. This explains why the relationship between 
interlace polynomials and Tutte-Martin polynomials discussed in [TJ [31 [T8] in- 
volves qjy and Q, but not q. The Tutte-Martin polynomials of isotropic systems 
cannot describe q because they involve arbitrary labels. 

The last polynomial we mention here is Courcelle's multivariate interlace 
polynomial [26]. If G is a looped graph with n vertices then C(G) is a polynomial 
in In + 2 independent indeterminates given by 



where GVB denotes the graph obtained from G by toggling loops at the vertices 
in B. Note that C(G) is essentially a table of the GF(2)-nullities of all the 
matrices obtained from adjacency matrices of full subgraphs of G by toggling 
some diagonal entries. Consequently C(G) contains the same information as 
the version of Q\(G) that uses the indeterminates in U {4> v , Xvi^Pv \ v £ 

V(F)}] as labels. This information is packaged differently in C(G), using two 
indeterminates and two possible loop statuses at each vertex rather than the 
three labels of Q\(G). The re-packaging obscures the basic theory of Q\ given 
in Section 3; [55] contains no analogue of Theorem [THl and the analogue of 
Theorem 1 171 is quite complicated. 




AUB\-u((GVB)[AUB]) v({GVB)[AUB]) 



Ans=0 



22 



7 Reduction formulas 



As discussed in [4j[9||30], evaluating q or q^ is #P-hard in general. Certainly 
the same holds for Q\, which evaluates to q and qw (with appropriate labels). 
In contrast, Courcelle [55] used techniques of monadic second-order logic 
to show that computing bounded portions of his multivariate interlace polyno- 
mial C is fixed-parameter tractable, with clique-width as the parameter. (The 
restriction to bounded portions of C is necessitated by the fact that 3' y ^ G ^' dif- 
ferent products of the indeterminates x a , y a appear in C(G).) The techniques of 
monadic second-order logic apply to a broad variety of graph polynomials, but 
they have the compensating disadvantage of producing algorithms with very 
large built-in constants. Consequently it is worth taking the time to investi- 
gate special properties of particular graph polynomials, which may be useful in 
simplifying computations. For instance, Blaser and Hoffmann |10j have used 
tree decompositions and GF(2)-nullity calculations to refine Courcelle's result 
regarding computation of bounded portions of C . 

Q\ and C determine each other term by term, so the results of Courcelle, 
Blaser and Hoffmann apply to Q\ too. In this section we prove a related result 
using Theorem [4TJ which gives formulas for the labeled interlace polynomials 
of graphs that possess split decompositions. These formulas extend results of 
Arratia, Bollobas and Sorkin [3 regarding interlace polynomials of substituted 
graphs. 

7.1 Pendant-twin reductions 

Ellis-Monaghan and Sarmiento [30 showed that the vertex-nullity interlace 
polynomial qw can be calculated in polynomial time for bipartite distance hered- 
itary graphs, i.e. graphs that can be completely described by two types of 
pendant- twin reductions [7]. Their argument involved the relationship between 
the Tutte polynomials of series-parallel graphs and circuit partitions of the 4- 
regular graphs that arise as medial graphs of series-parallel graphs imbedded 
in the plane [1HJ [301 1511 1M] • We extended this result to the two- variable inter- 
lace polynomial q and general distance hereditary graphs [64], by showing that 
a vertex-weighted version of q satisfies reduction formulas which allow for the 
consolidation of pendant or twin vertices into a single relabeled vertex; these 
reduction formulas are analogous to the series-parallel reductions of electrical 
circuit theory. Similar reduction formulas were also used by Blaser and Hoffman 
[9] in their analysis of the complexity of interlace polynomial computations. 

Qx also satisfies pendant-twin reduction formulas, which are of use in recur- 
sive calculations for distance hereditary graphs. 

Proposition 31 Suppose v ^ w are nonadjacent vertices of a labeled simple 
graph G, which have precisely the same neighbors. Then Q\(G) = Q\({G — w)'), 
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where (G — w)' is obtained from G — w by changing labels at v: 



4>(G-wy(v) = 4>g(v)4>g(w) + ^g(v)^g(w), 

X(g-w)>(v) = 4>g(v)xg(w) + Xg{v)<Pg(w) + xg(v)xg(w) ■ y 

+ Xg(v)^g(w) +iPg(v)xg(w), and 
^{G-w)>(v) = cj) G (v)ip G (w) + ^ G {v)<j) G {w). 

Proof. If G = I(F, C), the proof is indicated in FigureH) each configuration 
of v and w in F gives rise to a corresponding configuration in F — w. In general, 
we verify that Q\(G) — Q\((G—w)') by checking that each Gp matrix obtained 
by applying Definition [T2l to G has the same Gi 7 '(2)-nullity as the corresponding 
G'p, matrix. For the configurations involving in G this equality is obvious, as 
the two matrices are identical. For the other configurations the equality is not 
quite so obvious. For instance, if we add the first two rows of the first matrix 
displayed below to each row in the set containing Mn, we conclude that 



/ 1 1 \ 

Oil 

1 1 M n M12 

\ M 2 i M 22 / 



/ 1 1 \ 

Oil 

Mn M12 

\ M 2 i M 22 J 



Mn Mi 2 
M 2 i M 22 



This explains why ip G (v)ip G (w) is included in 4>( G ~w)'( v )- Similarly, 



/OOl \ 

1 

1 1 M n Mia 
\ Mai M 22 / 



/OOO ^ 

1 

1 Mn Mi 2 

\ M 2 i M 22 / 



1+v 






1 





1 


M n 


Mi 2 





Mai 


M 22 



explains why Xg(v)xg(w) ■ y is included in X(G-w)'( v )- ■ 

We leave it to the reader to verify the next two propositions, which give 
analogous results for adjacent twins and pendant vertices. For circle graphs, 
the propositions may be verified by completing analogues of Figure [8] for the 
two configurations shown in Figure |9] 

Proposition 32 Suppose v ^ w are adjacent vertices of a labeled simple graph 
G, which have precisely the same neighbors outside {v,w}. Then Q\(G) = 
Q\((G — w)'), where (G — w)' is obtained from G — w by changing labels at v: 

4>(G-w)<(v) = (j> G (v)(j) G {w) +Xg{v)xg(w), 
X(g-w)'{v) = 4>g(v)xg(w) +xg(v)<Pg(w), and 

i>{G-w)>{v) = 4> g (v)i/j g (w) + xg(v)^g(w) + -0g(w)0g(w) + iPg(v)xg(w) 
+ il) G {v)il) G {w) ■ y. 

Proposition 33 Suppose v ^ w are adjacent vertices of a labeled simple graph 
G, and w has no neighbor other than v. Then Q\(G) — Q\((G — w)'), where 
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Figure 8: One vertex replaces two vertices that give rise to nonadjacent twins 
in the interlacement graph. 



XX -X 

Figure 9: On the left, one vertex replaces two vertices that give rise to adjacent 
twins in the interlacement graph. On the right, one vertex replaces two vertices 
that give rise to a pendant vertex and its lone neighbor. 
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(G — w)' is obtained from G — w by changing labels at v: 

<f>(G-wy(v) = <Pg(v)<Pg(w) + XG{v)(j) G {w) ■ y + Xg(v)xg(w) 

+ Xg{v)iI> g {w) + ip G (v)(j) G (w), 
X(G~w)'{v) = <j) G (v)xG{w) +ip G (v)ip G (w), and 
^(G-w)>(v) = <I>g(v)iI> G (w) +i/j g (v)xc(w). 

For ease of reference we implement the removal of an isolated vertex w in a 
similar way, incorporating information about the labels of w in updated labels 
for a different vertex v. 

Proposition 34 Suppose w is an isolated vertex of a labeled simple graph G, 
i.e., w has no neighbor in G. Let v be any other vertex of G. Then Q\(G) = 
Q\((G — w)'), where (G — w)' is obtained from G — w by changing labels at v: 

<t>(G-w)'{v) = <Pg{v)Q\({w}), X{G-w)'(v) = Xg(v)Q\({w}), and 
^(G-wy(v) = ^g{v)Q\({w}). 

Suppose G can be reduced to a single vertex using Propositions [31] -[34] As 
noted in Corollary 5.3 of [3D] , such a reduction of G can be found in polynomial 
time, by searching repeatedly for isolated vertices, degree-one vertices and pairs 
of vertices v, w with the same neighbors outside {v,w}. In order to recursively 
describe the value of Q\ for an i?-labeled version of G, we apply the formulas 
of the appropriate proposition at each step. As each step involves removing a 
vertex, these formulas provide a description of Q\(G) in polynomial time. 

This description may or may not provide a polynomial time computation of 
Qx(G), depending on the computational personality of the ring R. For instance 
if R is a polynomial ring with three indeterminates for each v € V(G), then Def- 
inition[T4l includes contributions from 3'^°)' different products of monomials; it 
is impossible to explicitly compute such a large number of terms in polynomial 
time. For such a ring, Propositions - GH] are not really reductions in a prac- 
tical sense; they simply exchange combinatorial complexity (expressed in the 
structure of G) for algebraic complexity (expressed in the label formulas). On 
the other hand, if we are working over Z and using 0, x an d "0 labels that come 
from a small set of constants and indeterminates (like those used in obtaining 
q n , q or Q as instances of Q\), then we can determine Q\(G) by evaluating the 
indeterminates repeatedly in Z, and interpolating. Each individual evaluation 
involves only arithmetic in Z, which is computationally inexpensive, so for such 
a ring Propositions [3T] - [M] provide a genuine polynomial-time computation. 

7.2 Split reductions 

We discuss the following definition only briefly, and refer the reader to Cun- 
ningham [27] and Courcelle [211 US] [26] for thorough presentations. 

Definition 35 Let H and K be disjoint simple graphs, each with at least two 
vertices. Suppose S C V(H) and T C V(K). Then the join of H and K with 
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Figure 10: Every interlacement graph arising from the configuration on the left 
splits. 



respect to S and T is the graph (H, S) * (K, T) — (A, T) * (iJ, S) obtained from 
the union H U A by adding edges connecting all the elements of S to all the 
elements of T . The sets V(H) and V(K) constitute a split of (H,S) * (A, T). 
If H and A are labeled then the vertices of (H, S) * ( A, T) inherit labels directly 
from H and A . 

As an abuse of notation we will find it convenient to use * also when H or 
A has only one vertex. That is, for any simple graph G and any v G V(G), we 
may write G = ({v}, {v}) * (G — v, N(v)), even though {v} and V(G) — {v} do 
not constitute a split of G. 

Definition 36 Let G be a graph with a split G — (H,S) * (K,T). Then the 
split reduction of G with respect to H is the graph obtained from K by adjoining 
one new vertex h, with N{h) = T . That is, the split reduction of (H, S) * (K, T) 
with respect to H is {{h}, {h}) * (K,T). 

Note that if G = (H, S) * (A", T) and \ V{H)\ = 2 then either the two vertices 
of H are twins in G, or else one vertex of H is isolated or pendant in G. 
Propositions 13 1 1 — 1341 tell us that in each of these cases, the new vertex h of the 
split reduction may be labeled in such a way that the split reduction and G have 
the same Q\ polynomial. As we will see in Theorem 1411 the same is true for 
every split graph: (H, S) * (A, T) shares its Q\ polynomial with a reduced graph 
({hs}, {hs}) * (A, T), in which one appropriately labeled vertex hs replaces the 
ordered pair (H, S). 

To motivate this result, consider a connected 4-regular graph F that has 
a "4-valent subgraph" E. That is, there are precisely four edges connecting 
vertices of E to vertices of F — E; see the left-hand side of Figure [10] for an 
example. Let G be an Euler circuit of F, and let S (resp. T) contain every 
vertex inside E (resp. outside E) that is encountered exactly once on each 
passage of G through E (resp. F — E). Then every vertex of S is adjacent 
to every vertex of T in X(F, G), and these are the only adjacencies connecting 
vertices of E to vertices of F — E in I(F,C). That is, if H and A are the 
subgraphs of X(F, G) induced by V{E) and V(F) — V(E), respectively, then 
1{F,C) = {H,S)*{K,T). 
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Observe that in this special circuit partition P of F involves one of 

three possible "whole- E transitions" that reflect the connections in P involving 
the four edges connecting E to F — E. Comparing these to the connections in 
C, we obtain a "whole-l? transition label" of P with respect to C, correspond- 
ing to the sum of all the label products that represent choices of transitions 
at the vertices of E that are consistent with the "whole-i? transition" of P. 
These "whole- E transition labels" may be used to duplicate n(F, C) using cir- 
cuit partitions of the simplified graph obtained from F by replacing E with a 
single vertex. Equivalently, if we begin a computation of Q\(G) = n(F,C) by 
applying Theorem [T7] repeatedly to eliminate the vertices of E, then we can 
obtain the three "whole- E transition labels" by collecting terms. (This way to 
structure a computation - exhaust an appropriate kind of local substructure, 
and then collect terms before proceeding - was applied to calculations of knot 
polynomials by Conway 23 ; he called the 4-valent regions of knot diagrams 
tangles.) 

The following lemma of Balister, Bollobas, Cutler and Pebody will be useful. 

Lemma 37 (Lemma 2 of \Sp Let M be a symmetric matrix with entries in 
GF(2), and p a row vector. Then two of the three symmetric matrices 

( P ° tr m) > (o m) ' and ( P tr m) 

have the same GF '(2) -nullity, and the nullity of the remaining matrix is greater 
by 1. 

Although we will not require it, we might mention a sharper form of Lemma 
1371 proven in [69| : two of the three matrices actually have the same nullspace, 
and the nullspace of the remaining matrix contains the nullspace shared by the 
other two. 

Corollary 38 Suppose H is a labeled simple graph and S C V(H). Given a 
labeled partition P G V\(H), let Hp be the graph obtained from Hp by adjoining 
an unhoped vertex whose neighbors are the elements of SnV(Hp), and let Hp 1 
be the graph obtained from Hp by attaching a loop at the new vertex. Then two 
of the three numbers 

v{H§,),v(H P ),v{Hf,t) 
are equal, and the third is greater by 1. 

Proof. Let M be the adjacency matrix of Hp; then the second matrix of 
Lemma [571 has the same GF(2)-nullity as M. If p is the row vector whose ith 
entry is 1 if and only if the corresponding vertex of Hp is an element of S then 
the first and third matrices of Lemma [37] are the adjacency matrices of Hp and 
Hp e (respectively). ■ 

Definition 39 Suppose H is a labeled simple graph and S C V(H). The 
type of P G P\(H) (with respect to S) is 1, 2 or 3, according to which of 
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v{Hp),v{Hp),v{Hp l ) (respectively) is the largest. The labels of if (with re- 
spect to S) are the following: 

4>{h,s)= ( n f n xg(v)\ I n ^)W (hp) , 

Pev x (H) \ve4>(P) J \ve x (P) J \vei>(P) J 

of type 1 



X (h,s)= n m«) n n ^k^- 1 

Pev x (H) \ve<p(p) J \ve x (P) J \ve-4>(P) 

of type 2 



<p(h,s)= n m«) n n ^ p) - 

Pev x (H) \ve4>(P) J \ve x (P) J \vei>(P) J 

of type 3 

Observe that if we use <j>(H,S), x(H,S) and ip(H,S) as labels for a one- 
vertex graph {hs} then 

Qx({h s }) = 4>{H, S) + X (H, S)-y + tp(H, S) = Q X (H). 

Lemma 40 Suppose v is a vertex of K and either v £ T or the open neighbor- 
hood Nk{v) 7^ T — {v}. Then there exist a graph K' and a subset T" C V(K') 
such that \V(K')\ < \V(K)\ and ((H,S) * {K,T)) V X - v can be obtained from 
(if, S) * (K 1 , T") through some (possibly empty) sequence of labeled local comple- 
mentations. 

Proof. lfv£T then ((if, S) * (K, T)) v x -v = (H,S)* {K v x - v, T). 

If v is an element of T whose open neighborhood Nk(v) is not T — {v}, then 
((if, S)*{K,T)) v x -v = (if, S)*(K x '-v, (T-{v})AN K (v)). Here A denotes the 
symmetric difference and (if, S) denotes the graph obtained from if by toggling 
all adjacencies between vertices of S, and interchanging \ an d V' labels at every 
vertex of 5". As N K (v) j= T - {v}, there is some x G (T - {v})AN K (v). Then 



{(H, S) * {K\ -v,(T- {v})AN K (v))) x 

= (if, S) * {{Kl - v)l, ((T - {v})AN K (v))A(N K , (x) - {«})), 

because the labeled local complementation at a; restores the internal structure 
of if. ■ 

Theorem 41 Let if and if &e simple graphs with labels in R, and suppose 
S C V(H) and T C V(A"). TTien 

Q A ((ff, 5) * (ff, T)) = Q x (({h S }, {h s }) * (K, T)) 

where {hs} is the one-vertex graph with vertex labels cf>(H,S), x(ff, S) and 
i>(H,S). 
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That is, (H,S) * (K, T) may be reduced to {{h s },{h s }) * {K, T) without 
changing the value of Q x , so long as hs carries the appropriate label values. In 
the special case |V(if )| = 2, this reduction is one of the reductions discussed in 
subsection 7.1. 

Proof. If T = then (H, S) * (K, T) is the disjoint union of H and K, so 

Q X ((H, S) * (K, T)) = Q X {H) ■ Q X {K) 

= Qx{{hs}) ■ Qx{K) = Qx(({h s }, {h s }) * (K,T)). 

If V(K) = T ={v} then let (H, S) * (K, T) = G. By dehnition, Q\(G) is 

m- e f n ^(v)) ( n xGf(«) ) ( n 



+x(v) E n m*) n xg(«) n ^ \ y ^ 

pgvuh) \v£<P(P) J \vex(P) J \veip(P) 



-En n xa(v) n ^ \ v 

PeVx(H) W</>(P) / \ve X (P) J \vEip(P) 



= 4>{v) ■ {(j>(H, S) + X (H, S)-y + 1>(H, S)) 
+ X (v) ■ (tf>(H, S)-y + X (H, S) + ip(H, S)) 
+1>(v) ■ (4>(H, S) + X (H, S) + j>(H, S) ■ y) 

= Qx(({h s },{hs})*(K,T)). 

We proceed by induction on |V"(AT)| > 1, with T ^ 0. The argument is split 
into several cases. 

Case 1. If If has a connected component K' that does not meet T, then K' 
is also a connected component of both ({hs}, {hs}) * (K, T) and (H, S) * (K, T), 
so by induction 

Qx(({h s }, {h s }) * (K, T)) = Q x (({h s }, {h s }) * (K - K', T)) ■ Q X (K') 
= Q X ((H, S) * (K, T) - K') ■ Q X (K') = Q X ((H, S) * (K, T)). 

Case 2. Suppose every connected component of K meets T, and there is an 
edge vw in K with w T. We would like to apply the recursive step 

Qx{G) = Mv)-Qx(G-v)+^g(v)-Qx((GI)-v)+xg(v)-Qx(((G^)1)-v) (*) 

of Theorem Q7| to wandw, with G = (H,S) * (K,T). 
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T£v$T then G-v = (H, S) * (K — v,T), G\ — v = (H, S) * (K v x - v, T) and 
{{GDI) -v = (H,S)* ((-RT)a -v, t )- Th ese three equalities still hold if G is 
replaced by {{hs},{hs}) * (K,T) and (H,S) is replaced by ({hs},{hs}), and 
the inductive hypothesis applies in each case. We conclude that Q\((H,S) * 
(K,T)) = Q x (({h s h{h s })*(K,T)). 

If v S T the situation is more complicated, because local complementation 
at v changes the structure of H. However Lemma [40] assures us that each of 
the three values of Q A in © is of the form Qx((H, S) * (K', T)) with \V(K' 2 )\ < 
\V(K2)\, so once again we may cite the inductive hypothesis for each summand. 

Case 3. Suppose now that every connected component of K meets T and 
there is no edge vw in K with w £ T; then V(K) = T. If there is an edge vw in 
K then we use © again. This time though we require Lemma [40] only for the 
second term, because the two consecutive local complementations in {{G™) v x ) 
have no cumulative effect on the internal structure of H. 

Finally, if there is no edge in K then as V(K) = T, N(v) = S for every v S 
V(K). Consequently the vertices of K are nonadjacent twins in (H, S) * (K, T), 
and we can consolidate two of them into a single vertex using the formulas of 
Proposition [31] ■ 

The following definition will be helpful in discussing the recursive implemen- 
tation of Theorem I4T1 

Definition 42 The split width of a graph G, sw(G), is the largest integer that 
satisfies these conditions. 

1. sw{G) < \V(G)\. 

2. If G = (H,S) * (K,T) then sw(G) < max{\V (H)\ ,sw{{{h s },{h s }) * 
(K,T))}. 

We saw in subsection 7.1 that the reductions of Propositions I3T1 - l34l provide 
a recursive description of Q\ for graphs of split width < 2. In much the same 
way, Theorem [44] provides a recursive description of Q\ for graphs of split width 
< s, for each fixed value of the parameter s. The outline is simple. 

1. Given a graph G with sw(G) < s, find a suitable split G — (H, S) * (K, T) 
by searching for a subgraph H with |V(if)| < s, whose vertices fall into 
two subsets: V(H) — S (whose elements have no neighbors outside H) and 
S (whose elements all have the same neighbors outside H). The number 
of candidates for V{H) is polynomial in n = \V{G)\, because |V(-ff)| < s. 
For each candidate for V(H), there are no more than 2 s candidates for S. 

2. Use the formulas of Definitionl3l)lto determine the labels 4>(H, S), x{H, S), 
and ip{H, S). These calculations involve finding the GF(2)-nullities of no 
more than 3 S+1 different Gi 7 '(2)-matrices, with each matrix no larger than 
(s + 1) x (s + 1). 

3. Proceed to calculate Q\ for the reduced graph ({/is}, {hs}) * (K,T). 
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As in subsection 7.1, the computational complexity of this recursive de- 
scription depends on the nature of R. In a polynomial ring with three inde- 
pendent indeterminates for each vertex, the number of operations required to 
compute Q\(G) is clearly exponential in n, and the recursive description pro- 
vides a polynomial-time computation only for bounded portions of Q\(G). As 
sw(G) < s implies a bound on the clique-width of G (see Proposition 4.16 
of [2l]), this analysis is similar to Courcelle's result regarding computation of 
bounded portions of C for graphs of bounded clique- width [55] . 



In Z or Q, instead, arithmetic is computationally inexpensive, and we deduce 
the following theorem. Blascr and Hoffmann |10j have proven a similar result, 
regarding evaluation of C for graphs of bounded treewidth. 

Theorem 43 If G is a Q-labeled simple graph then the problem of evaluating 
Qx in Q is fixed parameter tractable, with split width as parameter. 

Polynomials like qw, q and Q, which are evaluations of Q\ in Z[x] or Z[x,y] 
rather than Q, can be determined by evaluating repeatedly in Q, and then 
interpolating. 
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